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Abstract 



Oi 
D . 

^ . The ratio of bulk to shear viscosity is expected to exhibit a different be- 

haviour in weakly and in strongly coupled systems. This can be expressed 
by its dependence on the squared sound velocity. In the high temperature 
QCD plasma at small running coupling, the viscosity ratio is uniquely de- 
termined by a quadratic dependence on the conformality measure, whereas 
in certain strongly coupled and nearly conformal theories this dependence is 
linear. Employing an effective kinetic theory of quasiparticle excitations with 
medium-modified dispersion relation, we analyze the ratio of bulk to shear 
viscosity of the gluon plasma. We show that in this approach, depending 
on the temperature, the viscosity ratio exhibits either of these dependencies 
found by means of weak coupling perturbative or strong coupling holographic 
techniques. The turning point between the two different dependencies is lo- 
cated around the maximum in the scaled interaction measure. 
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1. Introduction 



The observation that the quark-gluon plasma (QGP) formed in ultra- 
relativistic heavy-ion colhsions at the Relativistic Heavy Ion Colhder (RHIC) 
behaves almost like a perfect fluid [H, 0| , inspired numerous theoretical efforts 
to quantify the transport properties of the produced strongly interacting 
medium. In particular, various attempts have been proposed to calculate 
bulk and shear viscosities for QCD matter by means of different techniques 
(see [sj for a recent review). Likewise endeavours have been made to reliably 
extract information on the viscosity coefficients from RHIC measurements 
by continuously improving the phenomenological analysis of data (see [3] for 
a recent review). While currently heavy- ion collision phenomenology seems 
to focus on the shear viscosity as a relevant quantity, a firm knowledge of 
the bulk viscosity as well as other transport coefficients is equally needed to 
characterize the properties of the QGP concisely. 

The bulk viscosity C is a particular transport coefficient vanishing, within 
many models, in the non-relativistic and ultra-relativistic limits as well as 
for conformally invariant systems. In strongly interacting matter, the bulk 
viscosity to entropy density ratio C/s is expected to be small outside the 
deconfinement and chiral phase transition region. In the vicinity of the tran- 
sition region, however, it is expected to become larg e and even to diverge 
at a second order phase transition j^, 0, 0, H, [q], Iio| . Therefore, a precise 
knowledge of ( might be essential to understand and quantify the transport 
dynamics of the strongly interacting fluid in the critical region. 

The shear viscosity to entropy density ratio r]/s for QCD matter, on the 
other hand, is expected to be large, i.e. rj/s > 1, at large temperatures 
T, whereas in the vicinity of the deconfinement transition temperature 
it exhibits a minimum. Such a behaviour is known for a variety of liquids 



and gases [11|, Il2| and was conjectured [12|, Il3| to appear also in strongly 
interacting matter. 

The ratio rj/s was shown to exhibit a universal behaviour in gauge th eory 



plasmas at infinite 't Hooft coupling and infinite number of colours [ll|, Il4j . 
which led to the conjecture of a fundamental lower bound {rj/ s)kss > l/(47r) 
{Kovtun-Son-Starinets bound in natural units) for any physical system 



15| . A similar universal behaviour for the bulk viscosity or other transport 
coefficients is in contrast not known 1^, Nonetheless, in specific strongly 
coupled and nearly conformal theories that allow for a holographically dual 



supergravity description [18|, |19|, the ratio of bulk to shear viscosity was 
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found to behave as 

ch ~ - vi^ , (1) 

where the proportionahty constant is of order 0{1) and v1 = dP/de is the 
squared sound velocity expressed as derivative of the thermal pressure P 
with respect to the energy density e. Recently, a lower bound on this ratio, 
{C/v)b > 2 (1/A; — t>^) (Buchel bound), was conjectured for strongly coupled 
gauge theories in k spatial dimensions as being valid for all temperatures. 



where the holographically dual supergravity description is allowed [19[. This 



bound is exactly saturated as kinematical identity for all strongly coupled 



theories with holographic dual related to non-conformal branes [20 



In contrast, in scalar field theory [21 
with massive particles of a thermal fluid 
is uniquely determined by 



as well as for photons interacting 
22| . it was shown that the ratio (/rj 



C/v = l5{l-v!y . (2) 



^3 

In QCD, bulk and shear viscosities were calculated for large temperatures 



and small running coupling within kinetic theory [23|, |2J]. Under the re- 
laxation time approximation by assuming equal collision rates for bulk and 
shear viscosities [23| , these results give rise to a ratio (/r] that behaves para- 
metrically like in Eq. Thus, in perturbative QCD (pQCD) the bulk 

to shear viscosity ratio depends quadratically on the conformality measure 
Avg = (1/3 — fg). This implies that the Buchel bound is violated in high 
temperature QCD at small running coupling. 

Even though the precise holographic dual to QCD is currently unknown, 
one could expect from Eqs. (II]) and ([2]) that for deconfined strongly in- 
teracting matter the ratio (/r] undergoes a gradual change from the non- 
perturbative to the perturbative regime, rendering its dependence on the 
conformality measure from a linear to a quadratic one. It is therefore inter- 
esting to study the temperature dependence of the viscosity ratio from 
towards large T in order to understand how the transport properties of the 
medium are influenced by its thermal properties. 

In this work, considering the gluon plasma as composed of quasiparticle 
excitations with medium-modified dispersion relation, we analyze the proper- 
ties of the bulk to shear viscosity ratio. The equilibrium thermodynamics of 
such a quasiparticle model (QPM) was shown to describe successfully lattice 
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QCD results on the equation of state and related quantities for pure SUc{3) 



gauge theory j25|, |26|. To study the transport properties of the quasigluon 



plasma, we apply an effective kinetic theory approach [27[. From the results 
for bulk and shear viscosities obtained in this way, we show that in a quasi- 
particle picture the ratio (/r] exhibits indeed the above discussed behaviour 
with Avg-. At large temperatures, i.e. in the perturbative regime, the dom- 
inant dependence of (/r] on the conformality measure is quadratic, whereas 
near the deconfinement transition it renders into a linear behaviour. The 
turning point between the two is located in the vicinity of the maximum 
in the scaled interaction measure. Thus, on the level of transport proper- 
ties, a phenomenological quasiparticle model of the gluon plasma provides a 
systematic interpolation between the regimes of weak and strong coupling. 

2. Effective Kinetic Theory 

Assuming that the kinetics of quasigluons is appropriately described by 
an effective kinetic equation of Boltzmann-Vlasov type, the phase-space be- 
haviour of the single-particle distribution function / follows from (£ + V)/ = 
C[f], where £ is the Liouville operator, V is the Vlasov mean field term and 
C is the collision term. In local thermal equilibrium, the functional C = 0, 
and the thermodynamics of the QPM should be recovered jii], 2^. This is 



accomplished by / — )■ n{T) = d {exp{E'^ /T) — 1) ^, where n{T) is the Bose 
distribution function for gluon excitations with energy = 2 + n(T) 
and n(T) is the temperature dependent (but momentum independent) gluon 
self-energy. The degeneracy factor for colour and polarization degrees of free- 
dom is d = 16. Self-consistency of this approach dictates the form of the 
Vlasov term. In particular, this implies that in local thermal equilibrium V 
must be related to the temperature dependence of the self-energy in order to 



maintain thermodynamic self-consistency j29|, |30 



3. Ratio of Bulk to Shear Viscosity 

The viscosity coefficients can be derived from the above effective kinetic 
theory by assuming small deviations from local thermal equilibrium. The 
expressions for bulk and shear viscosities obtained within relaxation time 



approximation read |29|, l30|, |31 | 

d^p 
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C(T) 



1 
T 



(2vr) 



(4) 



where r denotes the relaxation time, which is related to the collision term C, 
anda = T\dn{T)/dT^). 

To quantify the transport coefficients for the quasigluon plasma from 
Eqs. Q and (j4]), we use for the gluon self-energy [25 



n(T) = -T^G\T) , 
where G'^(T) is the temperature dependent effective coupling, 



(5) 



11 In [A(T-T,)/T, 



i2 ' 



(6) 



defined such that at high temperatures the expressions for bulk thermody- 
namic quantities of SUc{3) gauge theory are recovered by the QPM. With the 
parameters Tg/Tc = 0.73 and A = 4.3. c f. |30l| . pure SUc{3) lattice gauge the- 
ory results for the equation of state |32[ |33[ are nicely described. The quality 
of this description is quantified by x^/d.o.f. = 1.4 ■ 10~^ for the lattice data 



from 32[ and by x^/d.o.f. = 1.5-10 ^ for the lattice data from 33|- Within 



this model, the leading terms of the squared sound velocity at asymptotically 



high temperatures, where 1 3> ^ \T{dG'^/dT) \ holds, read [23 



1 5 
3 ^48^ 



dT 



G\T)T 



dG\T) 
dT 



(7) 



which implies that reaches 1/3 only for asymptotically large T. 

From Eqs. ([3]) and (jlj) and by assuming an equal relaxation time for both 
transport coefficients, the ratio of bulk to shear viscosity in a quasiparticle 
description of the gluon plasma may be written in the form 



c 



15 {Av',Y-30Av![U{T)-a]v',iT) 



+m{T) - af {vl{T)f 



(8) 
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Figure 1: Scaled temperature dependence of the bulk to shear viscosity ratio Eq. ^ (solid 
curve) for 1.02 < T/Tc < 1.75 compared with available lattice QCD results from [s^ l 
(squares). The displaced error bars are determined from the statistical errors in ^ and 
ij |34| . The non-displaced error bars are obtained from the conservative upper and lower 
bounds on given in [3J] and the statistical errors in rj. 



with momentum integrals 

^'^^^ = I (l|"^^)[^ + d-MT)]^p'-' . (9) 

The momentum integrals receive dominant contributions from thermal quasi- 
particle momenta, \p\ ~ T, with proportionality constant of order 0{1). 
They follow, for all T > T^, the hierarchy > Xq > > 0. The 

terms in Eq. (|8]), which are proportional to [n(T) — a], contain only the 
temperature derivative of the effective coupling G'^(T), since the leading T- 
dependence cancels out in this combination. They would, therefore, not be 
present in the case of a temperature independent coupling g, i.e. in models 
with a quasiparticle mass M{T) with trivial temperature dependence of the 
form M{T) ~ Tg. Moreover, for constant M one finds E g . ([8 [) for the ratio 
C/r/, but with [n(T) - a] replaced by and Av^ = |3, Q. 

The ratio (/t] from Eq. is quantified in Fig. [T]as a function of T/Tc for 
a gluon plasma described by quasiparticle excitations with medium-modified 
dispersion relation and a momentum independent relaxation time r. For 
comparison, also shown in Fig. [1] is the bulk to shear viscosity ratio obtained 
from lattice QCD results in |3J] . The bulk and shear viscosities calculated in 



lattice QCD are fairly well reproduced within our quasiparticle model |30 



Noticing the large errors in Fig. [H the temperature dependence of the vis- 
cosity ratio in the quasigluon plasma is seen in Fig. [T]to be as well consistent 
with that found from the lattice gauge theory calculations. 

For T > 1.5 Tc, the viscosity ratio is entirely determined by the pQCD- 
like behaviour from Eq. (|2]). However, for T near Tc, the non-perturbative 
contributions in Eq. become significant, resulting in a reduction of the (/r) 
ratio. As evident from Fig. [T], for T > 1.02 the bulk viscosity is smaller 
than the shear viscosity in the quasigluon plasma. For T — >■ 1.02 Tc, the ratio 
(/•)] increases and reaches about 0.78. Quite similar values for this ratio were 



reported in holographic approaches [35|, l36|, l37|] and used in hydrodynamic 



simulations when studying the influence of a finite ( and t] on the elliptic 
flow fss']. Nonetheless, we note that in a different holographic approach 19 
a much larger value for C/r] in this temperature region was predicted. 

The viscosity ratio C/v Eq. ([8]) suggests apparently both, a quadratic 
behaviour oc (Af^)^, as found in the regime of small QCD running cou- 



pling 23|, and a linear behaviour oc Av^, which is in line with the dependence 



found in specific strongly coupled systems near conformality [la, [la, ll9| . A 
similar combination of dependencies was also pointed out in Glliil] as phe- 
nomenologically relevant relation. One should emphasize, however, that all 
other quantities entering the (/t] ratio in Eq. (|8]) depend on T = T(Av^). 
Therefore, Eq. (|8]) cannot be understood as a strict expansion series in powers 
of the conformality measure. Rather, for a quantification any temperature 
dependence has to be converted into a Af ^-dependence. 

The behaviour of T with Av^ can be obtained by inverting the relation 
between the squared sound velocity Vg{T) and the temperature. In the con- 
sidered quasigluon plasma, (T) follows as 



1 + 



TJijT) d (U{T) 



3 J2(T) dT \ T2 



(10) 



where the entering integrals read 



Ji 



dx- 



(ev^-l)Vv/^^ 2v/^^' 
2(x2 + z) (ev"^^ - l) 
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Figure 2: Left: Scaled temperature dependence of Av^ in the QPM. Right: The ratio (/rj 
from Eq. ^ (sohd curve) as a function of the conformaUty measure Av"^ . The dash-dotted 
curve shows a hnear tit (/t] — aAv"^ + (3 with a = 3.78 and (3 = —0.305. The dashed 
curve depicts the quadratic contribution in Eq. (|5]), 15 (Aw^)^. The dotted curve shows 
Buchel's conjectured lower bound for C/v- 



and 

(e^^'+^ - l) + z 

with z = n(T)/T^. For asymptotically large T, Eq. ( fTOl) simplifies to Eq. ([7]). 

The relation T = T(Af^) following from Eq. (|TOl) is shown in the left 
panel of Fig. El The right panel of this figure quantifies the resulting de- 
pendence of the ratio (/r] on the conformality measure Av^ in the interval 
0.01 < Af^ < 0.3. For small values of Af^, that is for large temperatures 
T > 1.5 Tc (see left panel of Fig. [2]), the viscosity ratio shows a quadratic 
behaviour with Af^. Indeed, for asymptotically large T, one finds with 
Eq. dl]), that At;2 ~ -5/{2ATr^)T\dGydT^), where \T^{dGydT^)\ < 1. 
Then, all terms in Eq. (|8]) are at leading order in proportional to 
(T'^dG'^ /dT"^)"^ , which implies in turn, that in this temperature regime all 
terms in Eq. (|8]) are at leading order proportional to (Af^)^. However, since 
T^l2/I-2 ^ ^ 1, the numerically dominant contribution to (/rj 

in Eq. ([H]) is given by 15 (Av^) , and Eq. (E]) is recovered. 

With increasing Av^, i.e. for decreasing T (see left panel of Fig. [2]), the 
bulk to shear viscosity ratio renders numerically into a linear rise with Av^ 
as seen in the right panel of Fig. [2j In fact, with decreasing T towards T^, 
the quadratic term 15 (Av^)^ in Eq. (|8]) is cancelled to a large degree of 
accuracy by the other terms, which in the regime 0.17 < Av"^ < 0.3 sum up 
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to approximately —15 (At>^)^ + ol At>^ + /3 with constant coefficients a. and /9. 
Thus, a hnear relation C/r] = a Af ^ + /3 emerges as evident from Fig. H] (right 
panel). The onset of this change in the behaviour with the conformality 
measure is located around T ~ 1.15 Tc. This lies in the direct proximity of 



the maximum in the scaled interaction measure [30 . 

Fig. 12] (rig ht panel) depicts, in addition, Buchel's lower bound on the 
ratio C,/^ [l9| (dotted curve). In our quasiparticle description of the gluon 
plasma this bound is satisfied also only for temperatures T < 1.15 Tc, when 
Afg > 0.17. We note, however, that this observation is based on the used 
approximation of a momentum independent relaxation time, which is the 
same for both viscosity coefficients. Therefore, it would be interesting to 
study how a momentum dependent r in Eq. (jS]), as proposed in f3l!|, would 
influence the quantitative comparison of Eq. ([8]) with the conjectured lower 
bound and the quantitative dependence of the viscosity ratio on At>f . 

This demonstrates that within an effective kinetic theory approach for 
quasigluon excitations, which is supplemented by information from equi- 
librium thermodynamics, both limiting behaviours of the bulk to shear 
viscosity ratio with Af^ known from the weak 23|] and strong coupling 



regimes |16l . Il8l . Il9l | can be found. Thus, the quasiparticle model for the gluon 



plasma provides a phenomenological method to describe non-perturbative ef- 
fects. Within such a model, one can effectively describe thermodynamics as 
well as transport properties of the gluon plasma in the whole temperature 
range from asymptotically large temperatures towards the critical region just 
above deconfinement. 



4. Conclusion 

Within our quasiparticle model for the gluon plasma, which is based on an 
effective kinetic theory, we calculated explicitly the bulk to shear viscosity ra- 
tio C/r/ in the relaxation time approximation. We have shown, that this ratio 
exhibits, depending on T, a quadratic or a linear behaviour with the confor- 
mality measure Af^ of the medium. Similar dependencies are known from 



small running coupling QCD [23[ and holographic [18|, Il9| approaches, respec- 
tively. While at asymptotically large T the quadratic dependence on Aw^ is 
dominant in our model, the linear dependence becomes more important when 
approaching Tc, governing effectively the qualitative behaviour of the viscos- 
ity ratio for 1.02 < T/Tc < 1.15. The turning point in the behaviour is 
located around the maximum in the scaled interaction measure. Therefore, 
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besides a proper description of equilibrium thermodynamics of pure SUc{3) 
lattice gauge theory, our quasiparticle model of the gluon plasma also de- 
scribes its viscosity coefficients, cf. [30|, and reproduces the dependencies of 



their ratio on Av^ as expected in the weak and strong coupling regimes, thus, 
providing a systematic link between the two. For momentum independent 
relaxation times, assumed to be common for bulk and shear viscosities, we 
find a significant dependence of (/rj on v^. This might imply consequences 



for viscous hydrodynamic simulations, cf. the studies in [38|, [391], in particular 



for LHC energies, where a larger temperature range is probed. However, for 
any phenomenological application one would have to include quark degrees 
of freedom into the model calculations. 
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